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This article provides a unified introduction to the theory of electron 
and proton spin polarisation in storage rings and it provides a common 
starting point for the written versions of the four talks that I gave at 
Monterey. 

1 Foreword 

Each of the four talks that I gave at Monterey had to do with spin polarisation in 
storage rings and accelerators and in each talk I covered the relevant and necessary 
aspects of the theory. Indeed, three of the talks essentially dealt only with theory 
and there was considerable repetition of the basics. If the written versions were to 
reflect the talks as I delivered them, there would again be repetition but there would 
also be an apparent lack of connection between those topics which were specific to 
each talk. Thus a reader who survived reading all four articles might still not have a 
solid view of the connections between the concepts covered. So it seems appropriate 
to provide a common introduction to the theory. That is the burden of this article. 
This also provides a suitable opportunity to present a synthesis of the various ways 
of describing the competition between polarisation build-up and depolarisation for 
electrons that I have come across or contributed to over the last decade. Moreover it is 
an opportunity to lay to rest some confusions that have crept into the subject. Owing 
to space limitations I will not attempt to maintain a high degree of mathematical 
rigour but aim instead to impart a feeling for the issues and for our current level of 
understanding. I shall refer to this article as Article I. 

The written versions of the talks themselves will be referred to as Articles II, III, 
IV, and V as follows: 

II Longitudinal electron spin polarisation at 27.5 GeV in HERA. 
( D.P. Barber for the HERA Polarisation Group ) 

III The permissible equilibrium polarisation distribution in a stored proton beam. 
( D.P. Barber, K. Heinemann, M. Vogt and G.H. Hoffstatter ) 

IV Unruh effect, spin polarisation and the Derbenev-Kondratenko 
formalism. ( D.P. Barber ) 

a Updated version of a contribution to the proceedings of the 15th ICFA Advanced Beam Dynamics 
Workshop: "Quantum Aspects of Beam Physics", Monterey, California, U.S.A., January 1998. Also 
in DESY Report 98-096, September 1998. 
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V The semiclassical FW transformation and the derivation of the Bloch equation 
for spin-1/2 polarised beams using Wigner functions. 
( K. Heinemann and D.P. Barber ) 



2 Introduction 

Spin behaviour in the electromagnetic guide fields of storage rings is dominated by 
two effects: 

• Spin precession 

• Spin flip due to synchrotron radiation emission []. 

In existing proton rings and those that will be built in the foreseeable future only spin 
precession is of significance since the synchrotron radiation power emitted by protons 
is negligible. However, as pointed out by Sokolov and Ternov in 1964, radiative spin 
flip can, for electrons, lead to a build up of polarisation IJ . This phenomenon is then 
commonly known as the Sokolov-Ternov (ST) effect []. At the time of writing, the 
only known practical way of obtaining a stored polarised proton beam is to inject 
a prepolarised beam provided by a suitable source and then accelerate it. Nev- 
ertheless another method has been suggested and I will comment on that in Article 
III. 

In the remainder of this article I will provide a unified overview of spin precession 
and spin flip and show how to arrive at an efficient description of their combined 
effect. 



3 Spin precession 

Spin precession for particles travelling in the electromagnetic fields in storage rings 
is most conveniently described in terms of the Thomas-Bargmann-Michel-Telegdi 
(T-BMT) equation § g §: 



4- S = ti AS 
dt 



where S is the 3-vector describing spin in the centre of mass frame and 
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The vector B is the magnetic field, E is the electric field and 7 is the Lorentz factor. 
The vectors f and f are the position and velocity and evolve according to the Lorentz 

b But it will become clear later that the distinction between the two can become blurred in storage 
rings. Indeed resonant spin flip in nuclear magnetic resonance experiments can be viewed either as 
flip due to photon absorption or precession by tt around an effective horizontal field. 
c In these articles statements made about electrons will also apply to positrons except for appropriate 
trivial sign changes in mathematical expressions. 
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equation. The quantity a = [g — 2)/2 is the gyromagnetic anomaly. For electrons 
a pa 0.0011596 and for protons a pa 1.7928. The other symbols used here and elsewhere 
have their usual meanings. The derivations of the T-BMT equation by its authors 
were purely classical in spirit. The derivation by BMT was based on the requirements 
of relativistic covariance. However, Thomas combined conventional notions of spin 
precession with the relativistic effect now called Thomas precession |3|, [J Q. Note 
that Eqs. (1) and (2) reduce smoothly to the usual nonrelativistic limit. To obtain a 
clearer view of the implications of the T-BMT equation one can rewrite it in terms 
of the field components perpendicular and parallel to the orbit: 

at mcy 

= — A ((f )S« + (- - 1 + § )B±) , (3) 
mc zj 7 2 

where for this part of the discussion the effect of electric fields has been ignored. 
Eq. (3) shows that for motion perpendicular to the field, the spin precesses around 
the field at a rate 1 + 07 faster than the corresponding rate of orbit deflection: 

8d S pin = (1 + a<y)58 orbi t = a^59 orbit + S9 orbit (4) 

in an obvious symbolic notation. This precession rate is strongly influenced by the 
Thomas precession. This is contained in the term I/7 — 1. For electrons (g ~ 2) the 
total precession is strongly suppressed. For protons (g ~ 5.58) the relative suppression 
is much weaker. 

However, 'spin' is a purely quantum mechanical concept. Moreover, we are not 
working in a regime where electron-positron creation and annihilation are important. 
Thus a two-component description of spin should suffice and one should therefore 
look for a Foldy-Wouthuysen transformation (Article V) of the Dirac Hamiltonian 
(containing a 'Pauli' term for the anomalous magnetic moment) appropriate for the 
semiclassical regime of a storage ring. By 'semiclassical' I mean that for the high 
energies involved it should only be necessary to keep terms up to first order in h. A 
Hamiltonian of the required type was already written down in 1973 by Derbenev and 
Kondratenko (DK) |J and takes the form | 

h dk _ h dk 1 q /r\ 

"op — "'op,orb ' 2 op °Z> ' \ / 



where: 



hon,orb ~ ° ^ op + ^op j (6) 



^Thomas also provided covariant forms for his equation. 

e The subscript 'op' is to remind the reader that we are dealing with operators. In this case they 
operate on two-component wavefunctions. The fields in h„p are external fields. The derivation of 
this Hamiltonian is the subject of Article V. 
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and: 

-j e mc 2 g — 2 -> ec(g — 2) 1 _ _ 

2mc J op 2 4m J op {Jop + mc 2 ) 

e fg — 2 mc 2 \ -> 

+ 2^U^ + J op (J op + mc 2 ) J (7r ° p A ^ + ' (?) 

and where 7r op and J op are defined as: 



7T p Pop — , Jop ^/c 7T p • 7T D p -|- Tfl C 



Thus the DK Hamiltonian consists of a purely orbital part of zeroth order in h 
and a spin part of first order in h. The orbital part resembles the familiar form 
of the classical relativistic Hamiltonian from the textbooks |7] and the spin part is 
reminiscent of a Stern-Gerlach (SG) dipole energy term. As will be noted in Article 
V, at second order in h~ this Hamiltonian gains just extra orbital terms. All in all, 
the DK Hamiltonian has a satisfying and physically transparent form. It is then no 
surprise that in first order in h~ the Heisenberg equation of motion (EOM) for the 
kinetic momentum n op is the Lorentz equation with an additional term for the SG 
force. It is also clear that in first order in h the Heisenberg EOM for the spin (h/2)a op 
is a precession equation with the same form as the T-BMT equation, Eqs. (1) and 
(2), since the operator Q op has a structure equivalent to that of Q in Eq. (2). In a 
wave packet approximation and at first order in h the Heisenberg EOM lead to the 
T-BMT equation for the expectation value (cr op ) (= the polarisation) and the EOM 
for the expectation value (vr p) of the kinetic momentum of a wave packet is again 
the Lorentz equation modified by a SG term ||. Thus we have now put the T-BMT 
equation on a firm quantum mechanical footing and have shown that it is the natural 
outcome of a semiclassical approximation. Moreover (see Article V), we know how to 
calculate beyond first order in h if necessary. Note that the magnetic SG terms differ 
from the familiar textbook forms for slowly moving particles but reduce to them at 
low energy: our terms contain Thomas precession contributions so that, for example, 
g/2 is replaced by g/2 — 1 + 1/7 = + 1/7. A detailed discussion on the SG terms 
in the DK Hamiltonian and on the SG forces allowed by covariance can be found in 
IP where the EOM are given a classical interpretation. See also Article III. 

The full Hamiltonian given by Derbenev and Kondratenko to include radiation 
effects is 

f-.dk _ idk 1 7 dk 1 1 dk (q\ 
lb tot ~ lb op ~r lb rad ' lb int \ a ) 

where hf^ d is the Hamiltonian of the free radiation field and where 

1) ft — * 

hfnt = e(<Prad Aad) + -(°op " &rad) (10) 

C Zi 

describes the particle-radiation interaction. The operator Q ra d has the same structure 
as rt, op except that the external field operators (denoted by the subscript 'op') are 
replaced with radiation field operators (denoted by the subscript 'rad'). 
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4 Spin distributions 



In the last section it became clear that to first order in h the centres of wave packets 
move (classically) according to the usual Lorentz force modified by a SG term and that 
the accompanying ((h/2)a op ) obeys the T-BMT equation. Thus for many purposes 
the particles and their spins can be treated as if they are classical objects and we 
are then in a position to move beyond single particles and to discuss classical spin 
and phase space distributions. Article V shows how to arrive at spin and particle 
distributions directly from the density operator. 

To construct a classical treatment one uses the correspondences: 

(fop) r , (if op ) -> 7T , (^<v) -> I* (11) 

where £ is a classical spin of length h/2. Then with the Hamiltonian: 

= (12) 

with 

^ b = J + e-4> (13) 
and the Poisson bracket relations Q: 

{r jt p k } = $jk , {r^rj = {p^pj = {r j} ^ k } = {P^Cfc} = , 

3 

= ^2 £ jkm £m , {j,k — 1,2,3) , (14) 

771=1 

and where semiclassically the r in Eq. (2) equals c 2 7r op /J, the Lorentz (modified by a 
SG term) and T-BMT equations emerge from the canonical equations of motion: 

r = {f,0*}, 7f = {vr, f) dfc } + — , £ = (15) 

Since storage rings and accelerators have accelerating cavities which subject the parti- 
cles to time dependent fields and since the magnet geometry is fixed, particle dynamics 
is best described in terms of the canonical coordinates u = (x, p x , z, p z , At, AE) where 
x, p x , z, p z describe transverse motion with respect to the curved periodic orbit and 
At, AE are the time delay relative to a synchronous particle (at the centre of the 
bunch) and the energy deviation from the energy of a synchronous particle respec- 
tively. The independent variable is now the distance around the ring, s. There is a 
corresponding (classical) Hamiltonian, correct up to first order in h, 

h = h orb + i-Q, (16) 

/If we were working to second or higher order in h we would use the Moyal algebra [jlO|, |ll| . In the 
present case of first order in H this simplifies to the Poisson algebra. 
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which enables the EOM to be written in canonical form and this is derived from f} 



<//,- 



by standard means fl^ ]. If the ring is distorted (see below), u describes the position 
with respect to the resulting closed orbit. 

We now make the idealisation that the beam phase space can be described in 
terms of a smooth continuous density, w(u;s), which is a scalar function of u and 
the azimuth s f\ It is normalised to unity. In the absence of dissipation and noise 
(e.g. due to synchrotron radiation) and ignoring the effect of the tiny SG forces on 
the orbital motion, w is constant along a phase space trajectory and obeys a relation 
of the Liouville type: 

dw ~ 

-7— = {h orb ,w\. (17) 

OS 

If the beam is stable, i.e. if w is the same from turn to turn, then it is periodic 
in s and we write it as w eq so that w eq (u; s) = w eq (u; s + C), where C is the ring 
circumference. 

Having assigned a phase space density to each point in phase space we now assign 
a polarisation P(u; s) to each point Q. P is the average over particles of the unit spins 
2£/H at (u; s). Since the T-BMT equation is linear in the spin and since in this picture 

the spins at (u; s) all see the same fl(u; s), P(u; s) obeys the T-BMT equation 

dP =? 

=n(u(s);s) AP . (18) 

as 

Because Eq. (18) describes precession, \P(u;s)\ is constant along a phase space tra- 
jectory. To make closer contact with the synchrobetatron motion, we can rewrite 



Eq. (18) as [y, |14 



dP ~ - ~ 

— = {h orb , P} + Q{u; s) A P (19) 
as 

which is analogous to Eq. (17) and assumes that P(u; s) is different iable in all direc- 
tions in phase space. Note that the polarisation of the whole beam as measured by a 
polarimeter at azimuth s is the average across phase space: 

PM-f*u *».)**.). (20) 

If the spin distribution is stable, i.e. if P(u; s) is the same from turn to turn, then 
P(u; s) not only obeys the T-BMT equation, but it is also periodic in s and we write 
it as P eq so that P eq (u; s) = P eq (u; s + C). We denote the unit vector along P eq (u; s) 
by h{u; s) []. This also obeys Eq. (18) and is periodic in s: h(u; s) = n(u; s + C). 
On the (periodic) closed orbit h(u; s) becomes n(0; s) and we denote it by n (s) Q. 

s Note that in Article V the phase space density is denoted by 'p'. 

h This is equivalent to associating a spin density matrix with each point in phase space. 
1 With respect to n(u; s) the spin density matrix at (u; s) is diagonal. 

J Many authors make no clear distinction between n and ho and many use the symbol n for uq. 
This can sometimes lead to confusion. In particular the original symbol for dn/dS (section 5.3) was 
"fdn/dj Q and some have erroneously understood "fdh/d^ to mean Eodho/dEo where Eq is the 
design energy. 
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Obviously n (s) obeys the periodicity condition n (s) = n (s + C), i.e. n (s) is the 
ring periodic solution of the T-BMT equation on the closed orbit. In general it is 
unique. 

In the foregoing I introduced the invariant (vector) spin field n(u; s) by appealing 
to physical intuition. The underlying assumption was that the field h(u; s), which is 
supposed to obey Eq. (18) over the whole of the beam phase space, not only exists 
but is smooth (to correspond with our expectations of the spin distribution in a real 
beam) and is unique. However, the situation is not quite so simple as I will now 
explain by describing some qualitative aspects of spin motion. 

If a circular accelerator only had vertical (dipole) fields, vertical spins would not be 
affected and n (s) would be vertical. Moreover, according to Eq. (3) a non- vertical 
spin would precess around n (s) 07 times per turn with respect to the (periodic) 
design orbit. I call the quantity 07 the 'naive spin tune'. It represents the natural 
spin precession frequency of this simple system. It increases by one unit for every 
rs 440 MeV (~ 523 MeV) increase in the energy of electrons (protons). But some 
rings contain vertical bend magnets so that the design orbit is not flat. The ring 
might also contain solenoidal fields of particle detectors. In these cases a periodic 
T-BMT solution, n (s), on the design orbit still exists but is no longer everywhere 
vertical and it is given by the real eigenvector (with unit eigenvalue) of the one turn 
(orthogonal) 3x3 spin transfer matrix for this design orbit f\. Indeed, for the HERA 
electron ring (Article II) no is made longitudinal at the east IP by means of spin 
rotators. The number of spin precessions around fio(s) per turn is extracted from the 



complex eigenvalues of the matrix |T9|, |2(J. We call this the 'real spin tune' or just 



the 'spin tune' and denote it by v svin . In general it deviates from aj [J. 

If the spin tune were an integer, the one turn matrix would be a unit matrix 
and n (s) would not be unique. This lack of uniqueness also manifests itself in 
extreme sensitivity to field errors. The quadrupoles and other magnets in storage 
rings normally have unavoidable small misalignements so that the periodic (closed) 
orbit deviates from the design orbit. Likewise the real n (s) deviates from the design 
no (s) since a spin on the closed orbit now 'sees' horizontal dipole components in the 
quadrupoles. There is also a small shift in the real spin tune. The angle between the 
two n (s)'s is roughly proportional to the amount of closed orbit distortion. But it 



becomes very large if the design spin tune is close to an integer [21, Q since the spin 
motion is then coherent with the imperfection fields. The spins are then said to be 
near an integer resonance (sometimes called an 'imperfection resonance'). 

Particle bunches in storage rings have nonzero transverse dimensions and energy 



^However, for a^O, the constraint n(u; s) — h(u; s + C) is obviously not equivalent to an analogous 
eigenproblem for n(u; s) since in general a spin at (u; s) set parallel to n(u; s) does not map into 
itself over one turn. Thus the naive algorithm based on a one turn map (e.g. see page 27 in |l5[ ) is 
incorrect; in general a 'n' constructed in that way would not obey the T-BMT equation everywhere 
along an orbit. As a result, a 'n' constructed in that way should not be used to obtain the vector 
dn/d5 needed, as in section 5.3, to describe radiative depolarisation of electrons (e.g. see page 52 in 
p6[). However, h can be obtained as an eigensolution of a modified eigenproblem [T^]. See also 
footnote |. 

'Actually, the complex eigenvalues only deliver the fractional part of the spin tune. The integer part 
must be found by following the spin motion for one turn. 
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spread and the motion of a spin, compared to that of a spin on the closed orbit, 

depends on the position in phase space via the u in Q(u; s). For particles circulating 
for many turns the total disturbance to a spin can grow to become very big if there 
is coherence between the natural spin motion and the oscillatory motion in the beam 
characterised by the spin-orbit resonance condition: 

v SV in = m + m x Q x + m z Q z + m s Q s (21) 

where the m's are integers and the Q's are respectively the horizontal, vertical and 
longitudinal tunes of the orbital oscillations. 

The integer resonances (\m x \ + \m z \ + \m s \ = in Eq. (21)) can normally be iden- 
tified with the imperfection resonances already mentioned and driven by the periodic 
imperfection fields along the closed orbit. We have absorbed their influence into a 
large deviation of n (s) from the design direction. The spin-synchrobetatron reso- 
nances (\m x \ + \m z \ + \m s \ ^ in Eq. (21)) (sometimes called 'intrinsic resonances') 
are driven by the quasiperiodic fields seen by particles executing quasiperiodic syn- 
chrobetatron oscillations about the closed orbit. The sum \m x \ + \ m z \ + \m s \ is called 
the order of the resonance. An imperfection resonance is then a zeroth order reso- 
nance. Although I have just been discussing the behaviour near resonance of arbitrary 
spins it should now be clear that n, which is a special solution of Eq. (18) constrained 
to be periodic, should, just like n Q (s), also show extreme behaviour near resonances. 
This is confirmed by the analytical structure and numerical output from the algo- 
rithms used for its construction [JHJ]. Near integer resonances in a distorted ring no 
deviates strongly from the nominal direction for the perfectly aligned ring and near 
intrinsic resonances the difference h(u; s) — n becomes large and increases with the 
synchrobetatron amplitude u and with aj 

For 27.5 GeV electrons in HERA (see Article II) the r.m.s. angle between n(u; s) 
and no (obtained by averaging across phase space) is just a few milliradians away 
from intrinsic resonances and about 100 milliradians very near such resonances. For 
protons at about 800 GeV in HERA (see Article III) on the '1-cr' torus this angle is 
typically 60 degrees unless Siberian Snakes are employed. Figure 1 depicts invariant 
spin fields h 'attached' to vertical betatron phase space ellipses for three different fixed 
energies but for the same invariant vertical emittance. Other examples are given in 
Article III. 

Although I introduced n via spin distributions, the history of n took a different 
course which provides more insight into its meaning and properties. It was first 
introduced by Derbenev and Kondratenko |25|, |6| in the process of obtaining action- 
angle variables for combined spin-orbit motion by 'diagonalising' the Hamiltonian in 
Eq. (12) and this aspect was further illuminated in f2B |. A similar approach can be 



used on the Hamiltonian in Eq. (16) fl2| , p7j , |2~6 , p8fl . I now give a rough outline of 
the basic ideas. 

It is assumed that the orbital motion is integrable and one makes an s dependent 
canonical transformation so that h or b is replaced by h or b = YjilnQili where the Qi 



m Note that the terminology 'intrinsic' and 'imperfection' must be used with care since synchrobe- 
tatron motion can also give rise to zeroth order resonance phenomena p3l E4| . 
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The effect of energy variation on the spin field 





— J 










Figure 1: A typical n-field at three energies, the second of which is close to resonance. 



are the three orbital tunes and the three Ii are the components of the orbital action 
vector I. Then one describes the spin motion with respect to a set of orthonormal 
axes (a 'dreibein') {hi(u\ s), n 2 (u; s), n 3 (u; s)) attached to each point in phase space 
and requires that £ • n 3 is a constant of motion. Clearly, n 3 (u; s) must be a solution 
of the T-BMT equation at (u; s) []. At the same time the dreibein is chosen so that 
hi{u] s) = hi{u; s + C) ( i = 1 — > 3) and so that £ precesses around h 3 at a constant 
rate relative to n\ and n 2 . The rate, denoted by u(I), should depend only on the 
actions I. The vectors n\ and h 2 are not solutions of the T-BMT equation. The 
vector n 3 has just the properties of the vector n introduced earlier. This choice of the 
dreibein, which amounts to a u and s-dependent rotation of the axes for describing 
the spin motion, is achieved by a suitably designed u and s-dependent canonical 
transformation which delivers a final Hamiltonian (correct to first order in K) with 



the 'diagonalised' form [26 



h = Sj lisQiii + l-nvl. 



spin 



(22) 



where I sp in = £ • ^3 is now an integral of motion, the spin action, and the Qi and Jj are 
the corresponding orbital tunes and actions. The Jj differ from the Ii by SG terms 
||26|| . Note that the concept of spin tune has now been generalised; instead of the 

closed orbit spin tune v spin we have a spin tune v{T) depending on the orbital actions 
(but not on I sp i n ) which differs slightly from v svin and which reduces to v svin for zero 
orbital actions 0. Now, in retrospect, the definition of resonance must be refined; 
we should really use i>(I) in Eq. (21) instead of ^ sp j n . It should now be clear why 
we sought a definition of spin precession rate, i.e. spin tune, which makes the latter 
independent of orbital phases and the azimuth s. Spin tune should tell us something 
about the degree of long term coherence between the spin motion and the orbital 
motion and allow us to express this coherence by means of resonance relations like 



"The angle between two T-BMT solutions following the same point on an orbit docs not change in 
time. See also footnote ^. 

°Note that for 7^0, v{I) cannot normally be obtained from a complex eigenvalue of the naive 
one-turn eigenproblem discussed in footnote wL 
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Eq. (21) (with u(I)). But if we work with a 'fake spin tune' such as that obtained 
from the one-turn eigenproblem (see footnote and p9|) which depends on orbital 
phases so that the 'fake spin phase advance' per turn varies from turn to turn, we 
can make no statements about long term coherence. With this redefinition of spin 



tune the dreibein (ni, n 2 , n 3 ) is unique except at spin-orbit resonances PBl BD| and 



by this uniqueness the vector n 3 is just the vector h introduced earlier except for a 
possible difference of sign! The exotic (unstable) behaviour of n near resonance is a 
manifestation of lack of uniqueness at resonance. 

Now I return to the questions of smoothness and existence of a n obeying Eq. 
(18). Since z> depends on orbital actions, n is potentially nonunique at almost all 
points in phase space because the resonance condition is satisfied almost everywhere 
if we include resonances of arbitrarily high order. Thus n might not be differentiable 
in all directions in phase space However, algorithms, both perturbative and 

nonperturbative, for constructing approximations to n are available [f2l|] (see also 



Article III ) and experience with calculating h by the author and colleagues seems to 
indicate that resonance effects rapidly become weak as the resonance order increases 
so that only a limited number of relatively low order resonances are likely to cause 
trouble. Therefore in the remainder of the article it will be assumed that the spin 
field n is a legitimate tool in practice. Nevertheless, the technical and interesting 
matters of existence and smoothness are under active investigation [R2J and knowledge 



gained from this study will be incorporated in our treatment of spin distributions. 
An example of the dependence of v on amplitude can be found in JHJ. 

An extension of the numerical work reported in |29| but carried out just before 
this article was completed indicates that v actually 'jumps over' resonant values as 



the orbit amplitude is changed fl33| , 34fl . That work is based on the 'stroboscopic 



averaging' algorithm in the computer code SPRINT [30|. But even more recent 



results from a new version of the SODOM algorithm |L7|, corroborates these 
findings. This implies, contrary to traditional expectations based on perturbation 
theory, that the spin-orbit resonance condition of Eq. (21) is never exactly satisfied 
in non-perturbative calculations. However, near to resonance, h still exhibits exotic 
behaviour. 

Although h and £ both obey the T-BMT equation they are very different objects; 
h is a function of the dynamical phase space variables but £ is a dynamical spin 
variable and by Eq. (14) the Poisson bracket {n,!;} vanishes. Now that we have a 
classical integral of motion for the spin, namely I sp i n , we recognize n as a phase space 
dependent semiclassical quantisation axis corresponding to the quantum observable 
(h/2)a op ■ n. We also see that the quantisation axis coincides with the direction of the 
equilibrium spin field. As we will see later (h/2)a op -n is a key object in the analytical 
theory of equilibrium electron polarisation and indeed it was originally introduced as 
an aid to calculating the electron polarisation |25|, ||. The analysis becomes more 
complicated if the orbital motion is nonlinear but in practice one tries to use an optic 
for which the nonlinear effects have been minimised and tries to restrict the beam to 
a phase space volume such that the motion is almost integrable. 

One last point on the virtues of n: a calculation of electron polarisation with the 
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computer program SODOM [17 which exploits h agrees well with a calculation using 
the Monte-Carlo spin tracking program SITROS which contains no notion of h. 

The material on spin distributions presented in this section is applicable both to 
electrons and protons. The application to protons is the topic of Article III so that 
for the remainder of this article I will focus on electrons and in particular on the 
modifications by synchrotron radiation to the concepts already presented. 



5 The effects of synchrotron radiation 

Synchrotron radiation (SR) emitted by stored electrons has three main effects: it 
determines the phase space distribution and it brings about spin polarisation owing 
to spin flip associated with synchrotron radiation (the ST effect) but the stochastic 
element of SR also causes depolarisation. Thus SR brings polarisation but it also takes 
it away! As we have seen already and as we will see below spin motion is irrevocably 
intertwined with the orbital effects. I will begin by summarising the orbital dynamics 
and then discuss the polarisation and depolarisation effects in detail. 



5.1 Orbital phase space 

Although SR spectra can be estimated by classical means || SR is fundamentally 
a quantum phenomenon; it consists of single photons so that one can only make 
reliable predictions by using quantised radiation theory. One then finds corrections 
to the classical spectrum [[HJ. The work of Huang and Ruth |S7j presented at this 
meeting is a good example of recent quantum calculations. 

Most of the SR in conventional storage rings is generated in the fields of the dipole 
magnets defining the design orbit. A quantum treatment for this case of the effects 



of SR on the orbital phase space distribution was carried out in 1975 [[38fl using the 
Hamiltonian of Eq. (10). I will return to this later but here I will follow another route 
which has the advantage of exhibiting the transparency needed for this article. 

Photon emission in the dipole fields is largely incoherent and detailed calculations 
show that one can consider the photons to be emitted over short distances of the 
order of p/7 where p is the orbit radius []. Furthermore in practical storage rings 
the energy loss per turn of a single particle is small compared to the nominal energy. 
Thus the dissipative effect is weak and for example in HERA (Article II) an electron 
at 27.5 GeV loses about 80 MeV per turn. Then for many purposes it suffices to 
describe the radiation reaction power p(s) from SR using a classical model in which 
smooth classical radiation reaction power p c i(s) is overlayed with a 'delta correlated' 
('white') stochastic component 5p(s): 

p(s) = p d (s) + 8p(s) , (8p(s)8p(s')) = R(Eq, K)5(s - s') (23) 

where the parameter R quantifies the intensity of the noise and depends on the design 
energy Eq and the curvature K 



p At this point I recommend the reader to consult the chart of time scales for electron dynamics 
in [^5). We will need this on several occasions. Indeed, an appreciation of these time scales is 
indispensible for understanding the physics of electron storage rings. 
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The equations for u of deterministic orbital motion derived from a Hamiltonian 
are then modified by inclusion of damping and stochastic terms and in the (very 
good) approximation that the photons are emitted parallel to the particle trajectory 
and neglecting interparticle interactions the resultant linearised stochastic differential 
equation describing motion with respect to the closed orbit can be used to construct 
the Fokker-Planck equation for the evolution of the phase space density JJ9], [4(], |4T], [42 1 
[| I write this as 

-£ = £ FP ,r b v> , (24) 
where the orbital Fokker-Planck operator can be decomposed into the form: 



^FP.orb ~ Gharri + Co + C\ + C2 ■ 

The term ChamU) is associated with the original symplectic (i.e. phase space 
density preserving) motion and it contains just first order derivatives with respect to 
the components Ui (i = 1...6). The operators Cq and C\ contain zeroth and first order 
derivatives and account for damping effects. The operator £2 contains second order 
derivatives originating in diffusion effects. 

A central property of Eq. (24) is that w(u; s) reaches equilibrium with w(u; s) = 
w(u; s + C) within a few damping times. At HERA at 27.5 GeV the longitudinal 
damping time is about 7 milliseconds ~ 350 turns ~ (design energy) / (energy loss 
per turn) Furthermore w(u;s) is gaussian and since the radiation effects are 

weak, w(u; s) is very close to being a solution of the radiationless transport equation 
Eq. (17) but with the radiation effects determining the beam size and causing a tiny 
ripple in the emittances [] as functions of s. 

Now that we understand the effects of SR on orbital phase space we can move on 
to spin. 



5.2 The Sokolov-Ternov effect 

Only a very small fraction of the radiated photons cause spin flip but for electron 
spins aligned along a uniform magnetic field, the ]{ and if flip rates differ and this 
leads to a build-up of spin polarisation antiparallel to the field. Positrons become 
polarised parallel to the field. The transition rates for electrons are 

5^ e 2 7 5 fr / 8 



Wn 16 ml&\pf \ + 5V3 
5v/3 e 2 7 5 fr A 8 \ 

w ^^WwK~^%)' (25) 

For positrons, interchange plus and minus signs here and elsewhere. 



'Restriction to linearised motion enables me to describe the chief qualitative features to be discerned 
without undue complication. 

r For electrons I define the emittance of a mode to be the r.m.s. action of the mode. 
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The equilibrium polarisation in a uniform magnetic field is independent of 7, 



St 



W n - W 

w n + w 



IT 



It 



5^ 



92.38% 



(26) 



For a beam with zero initial polarisation, the time dependence for build-up to equi- 
librium is 



where the build-up rate is 



P(t) = Pst [1 - exp H/t )] 



5\/3 e 2 7 5 ^ 



m 2 c 2 |p| 3 



(27) 



(28) 



The time Tq depends strongly on 7 and p but is typically minutes or hours. 

However, the fields in storage rings are far from uniform but since the system 
is semiclassical, Eq. (25), which was originally obtained from solutions of the Dirac 



equation, can be generalised and according to Baier and Katkov[44] for electron spins 
initially aligned along an arbitrary unit vector £ the transition rate is 



W = — 
2r 



1 



§<!■*> 



+ 



5^ 



(29) 



where s = direction of motion and b = (s As)/\s\. This is the magnetic field direction 
if the electric field vanishes and the motion is perpendicular to the magnetic field. 
The corresponding instantaneous rate of build-up of polarisation along £ is 



'bk 



T 
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(30) 



But instead of spin flip rates we really need an EOM for the polarisation itself and 
if we neglect the effect of stochastic (synchrotron radiation) photon emission on the 
orbit and imagine that all particles remain on the closed orbit (CO), the equation of 
motion for electron polarisation as given by Baier, Katkov and Strakhovenko (BKS) 

H HI 



is 



dP 



n co AP- 



P--s(P-s) 
9 V ' 



5V3 



b(s) 



(31) 



Note that the T-BMT term Q co A P appears here as the output of the radiation 
calculation itself. 

By noting that the characteristic time for polarisation build up is much larger 
than the circulation time 0, and integrating the BKS equation (Eq. (31)) one finds 
the generalised Sokolov-Ternov formula for the asymptotic electron polarisation in 
arbitrary magnetic fields along the closed orbit 



"Again, see Mj for a compilation of time scales. 
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3 § ds (Ms)-Ks))/\p(s)\ 3 ^ 

nks ~ 5Vs n ° fd8[i-i(Ms)-W]/\pW { } 

So the polarisation settles down aligned with n (s), the periodic solution to the 
T-BMT equation on the closed orbit. In rings containing dipole spin rotators (Article 
II) the polarization can usually not reach 92.38% since no(s) is then not parallel to 
the field everywhere. The corresponding polarisation build-up rate is 

Tbks ~ 8 mlc*cf dS |p( a )|3 • {66) 



The above formulae show that in the absence of stochastic motion the maximum 
attainable polarisation is 92.38% instead of 100%. Why should this be? 

At the simplest level the reason is clear: the probability for reverse spin flip is 
nonzero (Eq. (25)). Nevertheless, 'lay observers' often imagine that spin flip has 
something to do with spin's trying to reach the lowest energy level of the two levels 
of a magnetic dipole in a magnetic field and that once the spin is in its lowest level 
it will stay there. Then 100% polarisation would be achieved. Also, reverse flip by 
radiation emission would defy energy conservation. 

However, we are not dealing with spins at rest but with spins 'sitting' on relativistic 
electrons which already have quantised orbital energy levels so that the prohibition of 
reverse flip by energy conservation no longer applies. From Eq. (5) applied to electrons 
in a uniform vertical magnetic field it is clear that the energy change associated with 
spin reversal from up to down is (1 + 07 = 7(1/7— l+g/2)) larger than the separation 
of orbital energy levels hio c where u c is the angular frequency of the orbit. So one 
could naively imagine spin flip occuring without radiation but simply by a change of 
orbital energy level. A related phenomenon involving exchange of orbital and spin 



energy has been proposed by Derbenev [[47], [L2| while commenting on the possible use 
of transverse SG forces in storage rings. See Article III. 

Note also that the splitting of spin energy levels is not simply proportional to 
g/2 but contains a Thomas precession term, which indicates that the spin motion 
is coupled to the orbital motion. Furthermore, the average energy of a synchrotron 
radiation photon is tens of KeV. This is many orders of magnitude greater than the 
separation of spin levels. Moreover, photons emitted during spin flip tend to have 
higher energies than those emitted without spin flip. In addition, the polarisation 
does not reverse its sign with respect to the magnetic field at g — but at g ~ 1.2 
H fyj . This results from the fact that (I/7 — l + g/2) appears in the Hamiltonian 
hf*t (Eq. (10))instead of just g/2. 

Finally, it is interesting to note that the synchrotron radiation spectrum and the 



polarisation effects just depend on the curvature (i.e. the geometry) of the orbit [[H 



So the same effects could be obtained by using electric fields to bend the trajectory 
instead of magnetic fields. 

These comments should convince the reader that in the laboratory frame we are 
not dealing with a simple two level spin system. For further discussions relevant to 
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this topic the reader is directed to the articles by W. G. Unruh and J.D. Jackson in 
these proceedings and elsewhere . 



5.3 Radiative depolarisation 

The stochastic element of photon emission together with damping determines the 
equilibrium phase space density distribution. The same photon emission also imparts 

a stochastic element to the u in Cl(u; s) and then, via the T-BMT equation applied to 
spin motion in the (inhomogeneous) focusing fields and in a simple classical picture, 



spin diffusion (and thus depolarisation) can occur |48| . The polarisation level reached 
is the result of a balance between the Sokolov-Ternov effect and this radiative depo- 
larisation. In practice, depolarisation can be strong and it is therefore essential that 
it is well understood. 

But how can we calculate the equilibrium polarisation? After all, the polarisation 
at a point in phase space is the average of the unit spins 2£/h contained in a small 
packet of phase space at that point. Now, for protons, the phase space density is 
conserved along a trajectory so that no particles are lost from such a packet but 
for electrons their stochastic motion means that spins are continually diffusing from 
packet to packet. For the orbital motion one then employs a Fokker-Planck equation 
for the particle density. But polarisation is not a density so that it is not immediately 
clear how to proceed. Moreover the ST effect must be included so that an analogue 
of the BKS expression for stochastic orbits is needed. I will mention the best solution 
to this puzzle later but in the meantime I will follow a path which roughly reflects 
the way that estimates have been made in practice. 

A clue to the next step is contained in the above comments about the equilib- 
rium phase space distribution resulting from weak dissipation. There, the phase 
space distribution settles down to a distribution close to an invariant solution for the 
dissipationless problem but with the widths of the distribution determined by the 
radiation. Assuming that one has significant asymptotic polarisation the character- 
istic depolarisation time must be similar to the polarisation time, namely minutes or 
hours. Both are orders of magnitude larger than the orbital damping times. Thus the 
analogue for the polarisation would be that the direction of the equilibrium polarisa- 
tion at each point in phase space would settle down close to the equilibrium solution 
of the radiationless problem, namely h(u;s). Furthermore, the 'spin emittance' i.e. 
the average of I sp i n — £ ■ h at each point in phase space, would be independent of u 
and s. 

As has been customary I will now adopt these plausible notions as working as- 
sumptions that at equilibrium a) the polarisation is parallel to n(u; s) and b) the 
value of the polarisation is independent of u and s. In particular, it is assumed that 
the spin tune hardly varies across phase space so that there are no 'local' spin-orbit 
resonances and therefore no polarisation 'absorbers'. I will offer support for the first 
assumption at the end of this article but in the meantime some support for these 
assumptions comes from noting that by integrating the BKS equation along a deter- 
ministic synchrobetatron orbit the polarization settles down very nearly parallel to 
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n 



pE5| in analogy with the solution on the closed orbit (Eq. [32p [[ Furthermore, a 



study of a special but exactly solvable model of spin diffusion |TJj shows that far from 



resonance the polarization settles down asymptotically very nearly parallel to h. 

This picture was first proposed by Derbenev and Kondratenko ||. In the absence 
of radiation s n = (h/2)a op ■ n is conserved. But in the presence of radiation one has 

^ = - h [h d r k ad + ht t ,Sn}. (34) 

This is evaluated in the equations following Eq. (4.2) in || and by writing s = 
(h/2)a op the essence of the physics can be stated (very) symbolically in the form: 

ds n ds _ dn , . 

—^ = — -n + s-—. 35 
dt dt dt v ' 

The first term describes the rate of change of s n due to pure spin flip at a point 
in phase space (pure ST effect). The consequent build-up of polarisation is a 'spin 
damping' analogous to orbital damping. The second term describes the change in s n 
due to the fact that when a photon is emitted, the particle jumps without a change 
of spin to a new position in phase space where it finds a new n which will in general 
not be parallel to the h at the initial point. The projection of the spin on the n-axis 
has thus decreased stochastically so that s n diffuses in analogy with the diffusion of 
the orbital actions. This is where the depolarisation comes in. Thus the effect on the 
polarisation of the stochastic journey of a particle though phase space is accounted 
for by defining an appropriate quantisation axis at each point in phase space. Photon 
emission imparts both transverse and longitudinal recoils to the electron but since 
a photon is emitted typically within an angle I/7 with respect to the direction of 
the electron, the effect of the longitudinal recoil (i.e. the energy jump) dominates: 
the electron remains at almost the same point in x and z but can suffer a significant 
change in energy. Then by neglecting the effect of transverse recoil Derbenev and 
Kondratenko arrive at the following expression for the equilibrium polarisation along 
the axes h: 

$ ds h-(h-§)/\p(sW 

v dk = -— ^ (36) 



5v/3 fds((l-l(n.sf + %(§) 2 )/\p(s)[- 



where ( ) s denotes an average over phase space at azimuth s and 5 = AE/Eq is 
the fractional energy deviation from the design energy PL This formula differs from 
Eq. fl32| ) by the inclusion of the terms with dn/dS and use of n instead of Uq. The 
derivative dh/85 is a measure of the change of n caused by fractional energy jumps 
5 and its presence corresponds to the fact that the main consequence of a photon 
emission is a change in particle energy. The phase space average of the polarisation 
is 

P av ,dk{s) = Pdk (n)s (37) 

''very nearly' means that the angle between the polarization and n is much less than the angle 
between no and n. 

"This is sometimes written as #7/7 [||. See also footnote [J 
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and (n) s is very nearly aligned along n (s) (see the angle estimate below). The value 
of the phase space average, P a v,dk(s), is essentially independent of s. 

The effect of transverse recoil can also be included but contributes derivative terms 
(see Article IV, Eq. (2)) analogous to dh/d5 which are typically a factor 7 smaller 
than dn/dS and can usually be neglected [j50| |5lf . This point will be dealt with again 
in Article IV . 

In the presence of radiative depolarisation Eq. (^) becomes 



- 2 1 



"* " 8 m\,?C f \ \p(s)\> I 

Away from the spin-orbit resonances of Eq. (21) n(u;s) ~ n (s). But near res- 
onances h(u; s) deviates from ho(s) by typically several tens of milliradians at a few 
tens of GeV and the deviation increases with distance in phase space from the closed 
orbit. The spin orbit coupling function dn/d5, whose square (dn/d5) 2 in Eq. (36) 
quantifies the depolarisation, can then be large and the equilibrium polarisation can 
then be small. For example if \dh/d5\ is 1 the polarisation will not exceed about 57%. 

Note that even very close to resonances, \(n) s \ ~ 1: the phase space average 
polarisation measured by a polarimeter is mainly influenced by the value of Pdk in 
Eq. (37). 

The nice thing about this formulation is that a very complicated calculation of 
the effects of radiation has been distilled into a formula involving a few strange coef- 
ficients (emanating from the radiation theory) and a classical solution to the T-BMT 
equation, namely h whose behaviour encapsulates all of the depolarisation effects. 

To get high polarisation, one must have (dn/d5) 2 <C 1 in dipole magnets. If n is 
independent of the position in phase space, the derivative is zero: all points in phase 
space have the same quantisation axis and there is no depolarisation. But storage 
ring fields are inhomogeneous so that h varies across phase space. Thus the vector || 
depends on the optic of the machine. The optimisation of the optic required to make 
dn/85 small is called spin matching [PH| . This will be mentioned again in Article II. 
An example can be found in ||52|| . 

The term linear in dn/dS in Eq. ( |36"D is due to a correlation between the spin 



orientation and the radiation power |15|, [53|]. Alternatively, it can be considered to 



result from the interference between the two terms in Eq. (35). In rings where ho is 



horizontal due, say, to the presence of a solenoid Siberian Snake |]54 |, dh/d5 has a 
vertical component in the dipole fields. This can lead to a build-up of polarisation 
( 'kinetic polarisation) even though the pure Sokolov-Ternov effect vanishes. The rate 



1S T dk 



1 

1 



The expression for r^ k in Eq. (38) can be found from a purely classical calculation 
of spin diffusion by evaluating the effects of the second term in Eq. (35) or by other 
means PR [27]. I3I5L and indeed this was the first use for h [E| . Then we have a mixed 
calculation: the spin flip is described by quantum mechanics and the depolarisation is 
described by classical diffusion. But obviously kinetic polarisation will not be found 
by that route and the exotic resonance structure examined in Article IV would be 
missed. So it is clear that a quantum mechanical approach is needed to get the full 
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picture. An observation of kinetic polarisation would be a nice vindication of this 
viewpoint. 

The Derbenev-Kondratenko formula (Eq. (36)) has been rederived in a very ele- 



gant way by Mane |55[. He introduces the concept of generalised spin flip whereby 
he calculates the transition rates due to photon emission from 'spin up along 
h(x,p x , z,p z , At, AE; s)' to 'spin down along n(x,p x , z,p z , At, AE — A ph ; s)' where 
A p h is the energy of the emitted photon. One also needs the corresponding 'spin 
down' to 'spin up' rates. Then by requiring that the polarisation has the same value 
over all of phase space and imposing the constraint that the total generalised up- 
to-down rate equals the total down-to-up rate, and solving for the polarisation one 
arrives at Eq. (36)! In this formulation, the concept of depolarisation never arises! 
Instead one just has a statistical spin equilibrium. In a perturbative calculation of 
((dn/d5) 2 ) |55| one finds a series expansion in powers of emittances, i.e. in powers 
of h, in which each term contains a product of resonance denominators which impart 
the resonant structure to h. As the order of h increases along the series so does the 
order of the resonances contained in the successive terms. 

Apart from the calculation of Bell and Leinaas (Article IV), there are two further 
quantum calculations which should be mentioned, namely those of Hand and Skuja 



alt [5q| . They choose tiq as the spin quantisation axis. When a photon is emitted, the 
electron jumps to a new orbit. By writing the phase space coordinates as functions 
of the radiation fields and including damping phenomenologically, they calculate the 
rates of spin flip along n ]51| and obtain an expression for the equilibrium polarisation 
which is equivalent to the Derbenev-Kondratenko-Mane (DKM) expression in the 
limit in which the derivatives of n are evaluated in the linear approximation as in the 
SLIM formalism pi , |H| pp| . Thus they only find the first order spin-orbit resonances 



21 1 . Indeed, since their (quantum) representation of the orbit has a form similar to 



the classical representation in [J571] their final expression contains terms which are 



equivalent and identical in form to terms obtained in []57] in a model of classical 
spin diffusion. The fact that they only find first order resonances can then be traced 
to their use of just a first order perturbative calculation and the choice of h as 
quantisation axis. To find the higher order resonances one should use n. However, in 
a second calculation again using n , they calculated to higher order and although 
the outcome is not very transparent, the terms beyond the leading order in h contain 
high order resonant behaviour. But in the end the moral seems to be that it is more 
efficient to choose a quantisation axis (i.e. an unperturbed eigenstate), namely h, 
which reflects the physics and invest the numerical effort in working with this. In this 
respect the DKM formalism provides a clean practical framework in which to calculate 
higher order effects. Radiation fields can cause spin flip. Now we can see how, by 
treating the external fields experienced by the spin as functions of the radiation fields, 
spin precession can be regarded as spin flip as suggested at the beginning. 

We have seen that there are several ways to approach the estimation of the equilib- 
rium polarisation. In practice analytical calculations are carried out using the DKM 
formalism and this has been a successful and essential tool for predicting the main 
qualitative features of polarisation in electron storage rings. However, I now present 
an approach which is in many ways more satisfying. 
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5.4 Phase space and polarisation evolution equations 

Earlier, while discussing the difficulty of finding a Fokker-Planck treatment for the 
polarisation I promised further insights and they follow now. 

We have seen that the orbital phase space density w obeys the Fokker-Planck 
equation, Eq. (24). Then if the ST effect is ignored and it is recognised that spin is a 
passenger subject to the T-BMT equation it may be demonstrated using a classical 
picture |L3| [58J that the spin diffusion is described by the 'Bloch' equation: 

^ = £ FP , orb V + QAV, (39) 

where fl was defined by Eq. (16) and V is the 'polarisation density ' which is defined 
as 2/hx (density in phase space of spin angular momentum). The Bloch equation 
for the polarisation density is linear and it is universal in the sense that it does not 
contain the phase space density |14]]. It is also valid far from spin-orbit equilibrium. 



So the trick is to work with the polarisation density instead of the local polarisation 
P(u;s). It is simple to show that V(u;s) = P(u; s)w(u; s) and then by combining 
Eqs. (17) and (19) it is already clear that in the absence of radiation, V obeys the 
radiationless limit of Eq. (39), namely Eq. (40) below. The evolution equation for 
the local polarisation in the presence of radiation is obtained by combining Eq. (24) 
with Eq. (39). One finds that it has a more complicated form than Eq. (39) owing to 
the presence of the second derivatives in the £ 2 in £ FPorb - It a l so contains the phase 
space density so that it is not universal. But the local polarisation can always be 
obtained instead as V(u] s)/w(u; s). 

The ST effect can be included by adding in terms from the BKS equation multi- 
plied by the phase space density, together with some terms to represent the interfer- 



ence between ST effect and diffusion |L4| and in fact the full Bloch equation for the 
polarisation density and the Fokker-Planck equation for the phase space density can 
be obtained from quantum radiation theory |38[ . 
In the absence of radiation we obtain: 

-7T = {hor b ,V} + QAV . (40) 
OS 

From here it is easy to see, in analogy with the case of the phase space, that since 
the orbital damping, orbital diffusion and ST terms are very small compared to the 
remaining symplectic and T-BMT terms, the equilibrium (i.e. periodic) P(u; s) will 
indeed be almost parallel to h(u; s), at least away from resonances f\. 

The Bloch equation for the polarisation density is free from assumptions of the 
kind we needed earlier and in principle it allows us to calculate everything we need 
from scratch by looking at the beam as a whole instead of focusing on individual 
particles to begin with. It is clearly the best starting point for discussing radiative 



"The commentary in (5^1 on our calculations described in Article IV [ p0[ , contains the opinion that 
it is unsafe to employ fi as the appropriate quantisation axis. The above discussion should be an 
adequate response. 
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polarisation. Furthermore, the spin diffusion part (Eq. (39)) can be set up for any 
source of noise in the orbital motion — we just need the appropriate C FPorb . For 
example it could be applied to scattering of protons by gas molecules. 

Conclusion 

Spin polarisation in high energy storage rings is an exciting and exotic topic. I hope 
that the reader now has a solid overview of the status of our understanding and will 
pass on to Articles II, III, IV and V. An overview of the experimental aspects of 
electron polarization and plans for the future can be found in |60|, |6T| . 
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